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Abstract In this work exact expressions for the distribution function of the accumulated area swept 
by reunions and meanders of N vicious Brownian particles up to time T are derived. The results are 
expressed in terms of a generalised Airy distribution function, containing the Vandermonde determinant 
of the Airy roots. By mapping the problem to an Random Matrix Theory ensemble we are able to 
perform Monte Carlo simulations finding perfect agreement with the theoretical results. 
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1 Introduction 

Since the seminal work of de Gennes on simple models of fibrous structures and its popularisation by 
Fisher, the study of vicious walkers (namely, walkers that do not intersect) have attracted attention 
during the last two decades due to the wide range of applications in various branches of science and 

the statistical properties of the area swept by a Brownian excursion is a problem that was origi¬ 
nally studied by mathematicians. The Laplace transform of the distribution of the area (known as 
the Airy distribution function) was first computed by Darling and Louchard [71IT5]. The derivation 
of its moments together with the actual distribution were obtained by Takacs |29[I301[3T] . The Airy 
distribution function has appeared in a number of problem from different areas from graph theory and 
computer science as well as in physical systems modelling one dimensional fluctuating interfaces [M 
120] . applications in laser cooling nun] and fluctuations of sizes of ring polymers |21j . 

In this work, we focus on obtaining exact expressions for the accumulated area swept by N vicious 
Brownian motions. This paper is organised as follows: in section we define the process of vicious 
walkers and the stochastic quantities we are interested in. In particular, we show in section that 
the Laplace transform of PDF of the accumulated area can be expressed as a ratio of two N-particle 
propagators, which can be constructed as Slater determinants of single-particle states. This is a straight¬ 
forward consequence of using Quantum Mechanics formalism or, alternatively, of the Karlin-McGregor 
formula. The resulting expression is fairly general. In order to extract some properties we consider two 
type of processes: reunions and meanders. Moreover we consider two type of boundary conditions at 
a: = 0: absorbing boundary conditions (section® and reflecting boundary conditions (section]^. The 
inverse Laplace transform is resolved in sectionlQ and a formula for its negative moments is obtained 
in sectionFinally, in section]^ using a mapping between processes of vicious walkers to ensembles of 
Random Matrix Theory we are able to perform Monte Garlo simulations and compare with the exact 
formulas obtained in the previous sections. 
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2 Model definitions 

Consider the set of trajectories defined by N one-dimensional Brownian particles x{t) = 
with T e [0,T] such that 0 < xi{t) < ■ ■ ■ < xn{t) and let us denote as Mn[x{t)\ its probability mea¬ 
sure. Next, let us define as Ai the area swept by particle i: 


A,= 


i(r)dr. 


The set of areas A = {Ai, ..., Am) is a random fV-tuple whose joint PDF (jPDF) is given by viz. 

rx(T)=Xf N / rT 

lx{0)=xi 

Notice that the jPDF Q contains the PDF of the area swept for each curve i = 1,..., 

v%\a,t) = 


^x{T)=Xf ^ f 

Vn{A,T)= / 'Dx{T)MN[x{T)]T\d i / Xi{T)dT-A, 

Jx(0)=xi \Jo 


nOO 

■ / dAi ■ ■ ■ dAi_idAi^i ■ ■ ■ dAMV m{A, T) , 

Jo 


( 1 ) 


( 2 ) 


and, in particular, the PDFs of the area for the bottom and top curves denoted as Tn{A,T) = 
V^\a,T) and Bn{A,T) = 7^^^(A,T), respectively. It also contains the PDF of the total area, 
defined as 


N 


Sn{A,T) = 


dA, 


■dAMVN{A,T)^[A-Y,A^ 


(3) 


In this work we present exact formulas for the PDF Sn{A,T) for the accumulated area and compare 
the expressions with Monte Carlo simulations. A thorough analysis based on Monte Carlo simulations 
on the marginals (A, T) for i = I,..., N will be discussed somewhere else. 


3 Path Integral Approach 


Combining Q and @ yields 


.x(T)=Xf f ^ fT \ 

Sn{A,T) = / Vx{t)Jv[m[x{t)\5\A-^ / a;,(r)dT , 

Jx{o)=xi y "'O / 


with 


M[x{t)] = - exp 


N 


- E / 


dxi{r) 

dr 


CNlxir)] , 


where Z is the normalisation factor and Cm[x{t)] is an indicator function imposing the constraint of 
vicious particles 0 < a;i(T) < • • • < xm{t) on the positive part of the real line. Further, the Laplace 
transform of Sn{A,T), 


pOO 

SN{X,T) = {e-^^)s^= J dA SM{A,T)e-^^ 

has the following representation using path integral approach: 


5jv(A,r) 


rx(T)=xf 


/ 'Dx{T)MN[x{T)]exp 

-Xy / Xi{T)dT 

Jx{0)—Xi 

. i=l -^0 


(4) 
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In Quantum Mechanics formalism, the expression can be expressed as the ratio of two propagators 

G‘'^\xf,T\x,,0) 


5jv(A,r) = 


G^°\xf,T\x,,0) 


(5) 


with (a;, T|xo, 0) = (x\e '^\xq) and Hamiltonian operators 


N 


N 


N 


= E E = -b E 9 ". + E ’ “ = 0’ 1 


( 6 ) 


and a confining potential in M+, while = Xx for x > 0 and infinity for x < 0. Besides, 

Brownian particles are not allowed to cross which in QMf corresponds to having fermions. Considering 
the latter and that the potentials are one-particle type operators, the eigenfunctions of are given 
by the Slater determinant of the eigenfunctions corresponding to the one-particle problem /i(“^, that is 

(y, T|x, 0) = E , (7) 


with 




VM 


det 


N 

i=l 





Once the corresponding one-particle problems of the numerator and denominator have been solved, we 
have solved the problem via the formula ([^ . However, the corresponding expression is still too general 
to extract general mathematical properties. Thus, in light of the previous work we will particularise 
to the cases of reunions and meanders. 

For reunions we mean the set of processes by which the N vicious Brownian paths start and finish at 
the same point, that is Xi = Xf = 0. This is achieved by taking Xi = xj = e with e = (ei, - • • ,eAr), 
and taking the limit e to 0, viz. 


SP{X,T) = hm 


. GW(e,T|e,0) 


G^'(e,T|e,0) 


( 8 ) 


For meanders we mean the set of processes where the paths start at Xi = 0 and are allowed to finish 
at any final position Xf at time T. As such we need to integrate over all possible final positions and, 
therefore, the Laplace transform of the corresponding PDF reads 


rr\ r /wjv 

5b (A,r) = hm- J- - 


(9) 


where Wjv stands for the Weyl chamber Wjv = {x S M^|0 < Xi < • • • < x^} and d^x = dxi ■ ■ ■ dxjq- 


4 Absorbing Boundary conditions 

We note in the two cases above that we will need the one-particle problem for a = 1 (the one- 
particle problem for a = 0 is treated in the appendices). For absorbing boundary conditions (i.e. the 
corresponding eigenfunctions of the one-particle problem obey <j){x = 0) = 0) it has the well-known 
solution 




|Ai'(-a„)| 


Ai[(2A)i/3x- 


C^n] ) 


3 ( 1 ) = a„A2/32-i/3 , 


n = 1, 2,... 
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where Ai(a;) is the Airy function with zeros z = —with n S {1, 2,...}. Thus the eigenfunction for 
the A^-particle problem is: 




yM l<i,i<N 




( 2 A)^/® 




det [tfmizj)] 


( 10 ) 


where we have defined ipn{z) = Ai(z — an) with z = (2A)^/^a;. 


4.1 Reunions 


To obtain an exact expression for S^\x^T) given by (§ we need to perform the Taylor expansion of 
the Slater determinant (10) containing the Airy eigenfunctions. To do so we make use of the following 
formula m- 


^ oo ^ 

Ai(t)Bi(z + t)- Bi(t)Ai(z + t) = - V ^Qe{t), 

with Qniz) the Abramochkin polynomials m- In our case, we take that t = —an which yields 

°° z^ 

Ai(z - an) = Ai{-an) ^ -j^Qt{-an ), 


where we have used the identity 7 rBi(—a„)Ai^(—a„) = —1. Combining this together with the standard 
formula of the Taylor series of a Slater determinant (see appendix]^, we arrive at 


det ipr 


i{zj) = 


N 


n Ai'(-«",) 


b=i 


(2A) 


I 2^7 = 1 


te(7jv 




det [x, 


det [Qt- (-a„ 


:)], ( 11 ) 


where ^ S Un stands for the sum over all ordered fV-tuple indices ^ = (£i,..., ^yr), that is with 
< €2 < • • • < In- 

Everything boils down to be able to find an expression for the determinant deti<i(—)] of 
the Abramochkin polynomials evaluated at the Airy roots. It turns out that it is possible to obtain 
the lowest contribution in e if we recall that the odd-numbered polynomials are such that < 3271 - 1 - 1 ( 2 ) = 
z" -|- lower order terms, while even-numbered polynomials can be expressed in terms of the previous 
odd-numbered ones. This automatically implies that the first non-zero lowest-order contribution in e 
in the sum in eq. © is given by taking = 2{j — 1) -f 1 for j = 1,..., A^. This eventually yields: 




(2A)^ 


N 


Ai'(-a„ 


Li=i 


■ N 

n 


Xi 


An{xI, 


. . X 


N 


)Z\7v(an 


. CHrj 


,), ( 12 ) 


IV(iV —1) 

where we have defined ^n = (—1) ^ /[l! • • ■ (2A^ -I- 1)!], and Ajv(ai,..., ayr) is the Vandermonde 

determinant: 


AN^ai, ..., Cat) = 


1 • 

■ 1 

ai • 

• a-N 

N-l 

„N-1 

1 

■ Oat 


= n “ “j) ■ 


To obtain the expression (12 1 we have used that 

N 


det \x 


2(i-i)-sii 


.i-l 


N{N-1) 


det = (- 1 ) 2 /\^( 


l<i,j<N 


Otri\ t ’ • ’ 


)■ 
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The final expression of the Taylor expansion for the fV-particle eigenfunction is: 


iV(2iV + l) 




TT Xl{-an^) -pj- ^ ^^2 ,2 

11 |Ai'(-a„J| 11 * 


■ ■ ■ ^N^^Niptn-n ■ ■ ■ Oin^) I (13) 


with Pm = 'J n/V X fl. A similar derivation can be done for the denominator (see appendix [ b|. Using the 
results (13), ( [30| ), and ( [^ in the corresponding expression Q of the propagator, and after rearranging 
terms, we eventually obtain the Laplace transform of the PDF S^\A,T), viz. 


SI;\X,T) = ^ 




A' 


(r) 


E . 


aM _ 
— 


‘-AT ni,...,niv 

1 (.) JV-I 

2 2 ^JV -PT 


(14) 


rAf 


n A(2 + j)r(3/2 + j ), = l7V(2iV + 1). 


i=o 


4.2 Meanders 

For the case (|^ we report the final result (the derivation can be found in appendix [C|) 


g<rl(A.r) = 5: B'Pa„ .o„)4„(„„. 


Ei 


A 


(m) 

N ni,...,niv 


(15) 


where 


(o!ni, ■ ■ ■ , CHruv) jy ».// x ; 

11*=! Ai (-a„J JWj 


1 


[ det [Ai(a:^- - a„,)] , ^ . 

.Jwm i<uj<JV 3 


In this case we were unable to find a simple expression for the normalisation constant . 


5 Reflecting boundary conditions 

So far we have preoccupied ourselves by considering absorbing boundary conditions at a: = 0. We have 
also considered how the preceding derivations change when we take reflecting boundary conditions 
instead. Apart from the mathematical curiosity, the resulting process is interesting as it is related 
to the work m for processes ^/^’°)(t) which show a transition from type D to type D’ matrix 
ensembles (see [13] for details). 

The first thing is to notice how the one-particle wave-function changes in this case. If we denote as 
z = —Pn with n S 1,2 ,... the zeros of the derivative of the Airy function Ai (a:), then 


= //(!y?'| Ai((2A)''»» - M . P = /3„AV32-'/> 

Thus, in this case the A^-particle wavefunction takes the following form: 

( 2 A)^/® 


<F«(a;) = __ 


ViVinti v^|Ai(-/3„J| 


det [ipm [zj )], tpniz) = Ai(z - /3„), z = (2A)^/^a:. 
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5.1 Reunions 

We are left with deriving a Taylor expansion of the Slater determinant, as before. In this case, the 
other Abramochkin polynomials become useful: 


1 

{z)Ai{z + t) — A\ {z)Bi{z 1) = — —P^(z). 

TT n\ 

n—O 


From this expression we choose t = —fin to write 


Ai(z - Pn) = Ai(-/3„) ^ , 

where we have used the property 7rBi^(—,S„)Ai(—= 1. The expansion of the Slater determinant 
yields 


, det 


N 


nAi(-/3n,) 


.i=i 



det [x, 

l<i,j<N * 


. det ^(-/3 )], 


Again here, the lowest contribution in e is given by the even-labelled polynomials P 2 n{z) = z'^ + --- 
for n = 0,1, 2... and with Po{z) = 1. Thus we must take £j = 2{j — 1) for j = 1,..., A' yielding the 
following result for the Slater determinant 


tf'W(e)=5^(2A) 


N{2N-1) 

6 


TT 

i=i yS“|Ai(-/ 3 „,)| 


{X\ I ■ ■ • I (/3ni j ■ • ■ I /inn) 


N(N—1) ,r . . / s . /-Tn 

with 6]\f = (—1) 2 /[0!2! • • • 2(A— ^!v A!]. A similar analysis can be done to the propagator in the 

denominator of eq. ([^ (see appendix]^ eventually obtaining 




(Ar3/2) 




A 


(r) 


N 


■ ■ ■ j/Z^nn) /3„ 

/3ni • • • Pnn 


1 (r-) N-1 

(r) _ 22-^JV 


Ay) — 
— 


tN 


1=0 


+ A« = -A(2A-1), 


5.2 Meanders 

Similarly, in the case of meanders (see appendix]^, we obtain 


a».,, j,, _ (Arv>) 




A 


m) 

N 


o ^AJv(/3r^l,■•■/3r^«)„_(AT3/2)2/32-l/3y^«^^^ 
/ ^ \Pni^ • ‘ : PriN) R R ^ 

Pni ’ ■ ■ PriN 


( 16 ) 


where 


<)(/?„„... ,/3„„) = [ d^x det [Aiix, - /?„J] , = ^A(A - 1). 


and with no simple expression for the normalisation constant . 
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6 Inverse Laplace transform of .Sjv(A, T) 


First of all, we start by noticing that = Qisiis) with s = This, as pointed out already 

in [ 20 ], implies a scaling law of the form Sn{A,T) = so that QnIs) = £[QAr(a;)] 

with X = AT~^I'^. This applies to both cases of reunions and meanders. Notice that, even though this 
scaling appears mathematically, it can be easily derived by simple dimensional analysis. Indeed, as the 
dimensions of the diffusion constant are \D\ = L^T~^ and the dimensions of area in our problem are 
[yl] = LT (Nb. here our T refers to the time dimension, not to be confused with our final time T) then 
we have that 


SNiA,T) = 






as we have founc0 Secondly, by looking at the exponents and (see either at table or 

the expressions for the exponents in eqs. (14) and ( |16| )), it is clear that we must perform the inverse 
Laplace transform with S either an integer or one-third of an integer. 

To this end, we start from the result 


F{x)t 


^dx = 


„2/3 


F{x) = 2i/V(2i/3x-2/3)x-5/3 , 


with 

‘yij'iry 3 

J{x) = , 

and with U (a, b, x) the confluent hypergeometric function. After doing the rescalling of 
and a change of variables, we obtain the following formula 


2 - 1 / 3 ^ 52/3 


n+ 2 fc /3 -72 


or alternatively 


--v2-1/3,,2/3 


= / 

Jo 


nOC 

'd^+^F{x,'yy 

L 



e^^dx, F{x,^)=^J{^x , 




Qk-\-n 




(17) 


Here we have used the property 

n—1 

C[F^'^\x)] = s”F(s) - s^F^^-'^-'^\0 +), 




with notation F^'^\x) = dAF{x)/dx"', and the fact that any derivative of F{x,^) with respect to x at 
a; — 0 + is zero. 

To get an idea of the order of the derivatives involved in the function F{x,"f), one can see how the 
exponents and £1^^ vary with N and choose a value of the pair (n, k) which gives such an 
exponent. This choice is not necessarily unique and we show one possible choice in Table |l] The choice 
made is such that we only need derivatives with respect to 7 up to second order. The only case for 
which we are unable to use this prescription directly is for the case N = 1 for meanders, as this will 
imply the exponent = 1/3. A way around this is to consider in this case the pair (0, 2) which will 
give the derivative of the PDF, instead. 

This being settled, we proceed to find a simple expression for any order derivative of the function 
E(x,j). Starting from 

F{x,'y) = Q, ^,M(a;, 7 )^ , u{x,"i) = , 

1 For simplicity we have set D = |. This can be thought as equating dimensions of length squared with time, 
which implies that the dimensions of area are [A] = 
















# of Particles Absorbing Boundary Conditions Rejecting Boundary Conditions 


N 

p{r) 

n 

k 

Tp{m) 

-C//V 

n 

k 


n 

k 

p{m) 

n 

k 

1 

1 

1 

0 

1 

1 

- 

- 

1 

- 

- 

0 

0 

0 

2 

10 

2 

2 

0 

2 

2 

2 

0 

2 

0 

1 

3 

7 

7 

0 

3 

3 

0 

5 

5 

0 

2 

2 

0 

4 

12 

12 

0 

16 

A 

4 

2 


8 

2 

4 

4 

0 

5 


17 

2 

7 

2 

15 

15 

0 


6 

1 

6 

26 

26 

0 

12 

12 

0 

22 

22 

0 

10 

10 

0 

7 

35 

35 

0 

49 

15 

2 

91 

29 

2 

14 

14 

0 

8 

136 

44 

2 

Si 

20 

2 

40 

40 

0 


18 

1 

9 

ii 

57 

0 


27 

0 

51 

51 

0 

24 

24 

0 

10 

70 

70 

0 

100 

_3_ 

32 

2 

190 

_3_ 

62 

2 

30 

30 

0 


Table 1 Exponents and as a function of N for both absorbing and reflecting boundary conditions. 

In the columns n and k, we show the number of derivatives involved for the function E(x, 7 ). Recall that the 
expressions for the exponents are = N{2N + l)/3 and /3 for ABCs and E^^ = N{2N — l)/3 

and E^' = N{N - l)/3 for RBCs. 


and using properties of the hypergeometric confluent function we arrive at the following result (see 
appendix 


a'=+"E(a;,7) ^3 


n k y . X 

--^-s,-,u(x,7) , (18) 

VO o / 




1=0 s =0 


where the set of coefficients and are given by 


CT> = 


n\ 


2"-2^(n-£)!(2£-n)! ’ 


(£) = 1, (£) = -^(l + 2£) , (£) = -3 

dI,^\£) = ^(1 + 2£)(5 + 6£), d[^\£) = 3(7 + 6£), D^^>(£) = 9. 


' = 0,...,n, 


(1)/ 


1(2)/ 


Notice that for the set of coefficients {£>i^\£)} we have already taken into account the fact that we 
only need derivatives with respect to 7 up to second order. 

With the help of eqs. 0 and ( |18[ ), we can now perform the inverse Laplace transform of Q^\s) and 
Q^^\s), obtaining the following generalised Airy distributions for absorbing boundary conditions: 


N 


Qn\x)= , 

/ N 

Q^n\^)= X] S^Vni>---anN)^Ar(ani:---aniv)2^^’(™) 


ni. ,niv 

For reflecting boundary conditions we have instead 


Qn\x) = 


• Pun) q-{n,k) 


N 


...R 

ii Hun 


i=l 




>7)= y: B<,” ^V ffi, (x,f /i„ 

Pni ’' ' Pun V • 1 

ni,...,nAr ' 1 ' jv Y 


(19) 

( 20 ) 

( 21 ) 

( 22 ) 


In both cases we have defined 
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with s S {r{eunion),m{eanders)}. In the expression of 7 ) the pair of indices (n,k) must be 

chosen according to the values related to the exponents and appearing in table Notice, in 
particular, that for TV = 1 and (n, k) = (1, 0) we have that An = 1, = 0, = 2. 

Thus we recover to the so-called Airy distribution, viz 


n—1 





2/3 


2^3 

e 27x2 JJ 


54 273 \ 
6’3’27x^J ■ 


7 Moments 


In this section we discuss the derivation of the moments. We consider negative powered-moments. Let 
us start by fixing some notation. Let us denote the n-th moment for a PDF F(x) with x G 
An = /g°° dxx^F[x). Suppose next that the Laplace transform of a function F{x) is F{s) = s^e~ 
with E and a two constants. Then one can show that (see appendix [f| 


, 2/3 


= ^ ^ + ^)) 

2 F{:.) ■ 

In our case let us denote the moments as dxQ^^\x)x"' with a G {r, m}. Using the previous 

results, we obtain the following formulas for the negative moments of the PDF for reunions and 
meanders with ABCs 


nxF) _ 9- 

-'"jV.-i/ — ^ 


_ 9 - 




■E 




-1 


3r(liu + E 


N 




ni. ,nN 




Z\7v(ani, 


ni 


1 



,)Bn{c 





— I ) 


Similar expressions can be found for RBCs. 


8 Monte Carlo Simulations 


To check the correctness of our analytical hndings, we have performed Monte Carlo simulations exploit¬ 
ing the connection with Random Matrix Theory mm to generate samples of non-colliding paths. 
Following the notation in m we first recall the definition of the Pauli matrices 


tJi = 





while we denote the 2x2 identity matrix I 2 as cto- For meanders we define the following 2N x 2N 
Hermitian matrices (t) and (t) corresponding to absorbing and reflecting boundary conditions 

(0 = ^ (^) O) C) CTo -I- (t) ® ai+ (t, O) (g) (72 + (t) (g) CT 3 , 

“T (0 = ^ (^> O) C) CTo -I- JOy ^ (t, O) (g CTI -b (t) (g 0-2 -b {t) g (73 , 

respectively. Here the notation is a bit involved but it means the following: (•••(tH, O) stands for 
bridges starting at the origin and hnishing at the origin at time T, while and s^°‘> denote N x N 
antisymmetric and symmetric matrices, respectively, whose elements are either bridges or standard 
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Bottom Curve 






Fig. 1 Two instances corresponding to a reunion (top panel) and a meander (bottom panel) of N Vicious 
paths with reflective boundary conditions at a; = 0, generated by the method explained in the text. We also 
show the corresponding areas for the bottom, one of the inners, and the top curves. 


Brownian motions. Being more precise, if we denote as b(t) a standard Brownian motion with 6(0) = 0, 
and with B(t) a Brownian bridge with B(0) = B(T) = 0 then 


and 


(a) 


(t) = 


b^fit)/V2 i < j 

bu\t) i=j , 

< i 


(a) 


it) 


blf{t)/V2 i<j 
0 i = j 

-bji\t)/V2j < i 


Srljit^O) = < 


Blj\t)lV2 i < J 

Bu\t) i=3 


(a) 


{t,0) 


B\f{t)/V2 i<j 
0 i = j- 

-B^;\t)/V2j<i 


Here the index a = 0,1,2,3, simply states that we need to construct different matrices. Similarly, for 
the case of reunions we have instead the following two 27V x 2N Hermitian matrices for absorbing and 
reflecting boundary conditions 

“T it) = it, O) 0 CTo + Sy ^(t, O) 0 (Ti + Sy^(t, O) 0 (72 + {t, O) as , 

“T it) = it, O) 0 (7o + ia)p {t, O) 0 (Ti + foy ^ (t, O) 0 (72 + s® (t, O) 0 (73 , 


respectively. 

Bridges are easily generated from a standard Brownian motion. In its discrete version if R(k) is a 
standard random walk with k G {0,... ,iV} with i?(fc = 0) = 0 then a discrete Brownian Bridge is 
given by B(k) = R{k) — {k/K)R{k). 

To estimate the different PDFs we have taken K = 10^ and we have subsequently generated 
Af = 10® S{t) matrices for the four different process. For each instance s = 1,..., A/” the area swept 
by each path is estimated as 


i = l,..-,N. (23) 

k=0 

where A(fc) = (—X^ik), ..., —Ai(fc), Ai(fc),..., Xffik)) are the eigenvalues of the matrix S. This area 
is then rescaled to the cc-variable j. Finally, the sample set for each area 
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Fig. 2 Two instances corresponding to a reunion (top panel) and a meander (bottom panel) of N Vicious 
paths with absorbing boundary conditions at ® = 0 and generated by the method explained in the text. We 
also show the corresponding areas for the bottom, one of the inners, and the top curves. 


i = 1,..., N, is used to estimate the PDFs as: 


:,ii) 

N 


J\f 


(24) 


S=1 


Here we focus on the PDF Qn{x) of the accumulated area. An instance of the four processes can be 
found in figures [l] and 

Results of the PDFs estimated by Monte Carlo simulations and comparison with the theoretical for¬ 
mulas for the PDF Qn{x) of the accumulated area are reported in figure]^ for both types of boundary 
conditions. 


Looking at these figures some points are in order: the comparison between the set of formulas (191 


and (211 and Monte Carlo estimates could be done up to iV = 10 and N = 9, respectively, as otherwise 


the numerical evaluation of these exact formulas becomes prohibitively long. For the formulas (201 


and (22) we have the an additional numerical problem due to the evaluation of the and, as a 


consequence, the comparison with Monte Carlo simulations is performed up to iV = 2 in both cases. 


9 Conclusions 

In this work we have generalised the Airy distribution function of the area swept by one Brownian 
particle performing a reunion to the case of the accumulated area swept by N vicious Brownian 
particles. The exact formulas have been contrasted with Monte Carlo simulations showing perfect 
agreement. 

There are several open problems which we are currently investigating. First of all, whereas obtaining 
exact expressions for negative moments of the distribution of the accumulated area swept by N vicious 
walkers is rather simple, we wonder whether it is possible to obtain exacts expressions for the positive 
moments, perhaps generalising the techniques used for the case of one walker that can be found, 
for instance, in [2^I^[T4] . Secondly, we ponder whether it is possible to find exact formulas for the 
distribution of the area swept by either the bottom or the top paths. Within QMf this entails to being 
able to construct A^-particle eigenfunctions for distinguishable fermions. As this seems to be a daunting 
task, at the very least it would be interesting to study the properties of these distributions numerically. 
This can be done fairly efficiently, in particular for rate events, by combining the mapping to RMT we 
have used here together with the Wang-Landau algorithm, a numerical technique that has been used 
in a similar context to study the statistics of extreme eigenvalues [25j . This analysis can obviously be 
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q'-n\x) 






Fig. 3 Plots for for N = 1,2,3,4 and 10. The top plots corresponds to absorbing boundary condi¬ 

tions, while the bottom plots are for reflecting boundary conditions. Thick solid lines the theoretical results 


given by (19l. All PDFs have been estimated by generating 10“ samples. 


extended to the cases in which the distribution of jumps of the Brownian motion comes from either 
thick or thin tails. 
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A Taylor expansion of Slater’s determinant 

In this secion we discuss the Taylor expansion of a Slater determinant. Starting form the definition of the 
determinant and doing a Taylor expansion we write: 


^ sign(P)]^(p„^(i)(2:i) 
P^Sn ^=1 


lx In ^ 

Y. “en(p) E xTAxtnpi'Aw 


PGSn 


= E 

ii,-.- ,£n ^0 


£l! • • • £n'- l<i,j<N ^ 1 ^ ' 


(25) 


Due to the antisymmetric properties of the determinant, only those terms in the multiple sum with different 
values of £’s are different from zero. With this in mind we rewrite the sum over the ^ = (£i,... ,^iv) as a sum 
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over over Un of ordered indices (let’s say li < £2 <■■■ < In) and then we sum over its permutation, viz. 

An 




£€ Un 


= E 

IGUn 



£1 

^1 •• 

?JV>0 

4!-- 

1 

4!- 

•4jv! 

1 

4!- 

■■£nI 

1 

4!- 

■ ■ £n^. 


P^Sn 


PGSn 


(26) 


det A] det (0)], 


as we wanted to show. 


B Reunion 


For the normalisation factor we need to derive the propagator of N vicious walkers in the semi-infinite line. 
This propagator is this case can be written as follows 

G.(», Tk, 0) = jTj (I) " J,, -. - ^i) V4-’ , (27) 

where = (l/-\/]vT) deti<ij<iv[sin(giXj)]. Using the Taylor expansion formula for the Slater determinant 

we arrive at 




det det kf ■ 

- • Tvr <■ t J 1 ^ Tvr ‘ ' ' J 


VM (2si -I- 1)! • • • {2sn + 1)! i<i,i<iv^ ‘ i<i,j<N 
with s = (si,..., Sn). Noticing that the lowest contribution is obtained by of Si = i — 1, we obtain 


/ \ 2(^1 




AN{ql,...qN), 


(28) 


(29) 


iV(Ar_i) _ _ ^—_ 

with Pn = (—1) 2 /[l!. •. (27V — l)!v7V!]. At this stage we notice that the propagator goes like 


Go(e,r|a3,e)~/3^r-5^‘'‘^+^^ 


H' 






B, 


Bn = 


POO PO 

Jo Jo 


dqi ■ ■ ■ dqNC 




(30) 


‘^iv(9ii • ■ ■ : 9Ar) 


where the constant Bn is derived in appendix [P] 

For reflecting boundary conditions the expression for the propagator is the same but with the Slater determinant 
given by = (l/-^/)^) deti<ij<jv[cos(qiXj)]. Using the expansion formula for the Slater determinant we 

arrive at 




(2«i)! - (2™)! ■susd’’’ '• 


with s = (si,..., Sn). Noticing that the lowest contribution is obtained by of Si = i — 1, we obtain 


(^) ^PnT-^^^An {xl,...,xl) 




Z\]v(g?,...gw), 


(31) 


(32) 


iV(Ar_i) _ _ _ ^—_ 

with Pn = (—1) 2 /[O! • • • (2(7V — l))!v7V!]. At this stage we notice that the propagator goes like 

„ iV( 2 iV-l) 0/0 o \ 

Go{€.^T\x^e) Cal j5NT 2 yei,... Fn , 

/ n \ -^ POO PO 

I -L 

where the constant Fn is derived in appendix [D] 


(33) 


dqi---dqNe AN^qi,-■ ■ ,qN) 
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C Meander 


For the case of a meander of N vicious walkers (a star with a wall) we need to integrate over the final position 
X. Using the Taylor expansion of the Slater determinant as explained in the main text we can write 


f d^xG^j^\x,T\€,0) = f d^x4>^^\x)4>^^\e)e 

^ JW^ 


~ (2A) 3 yjY 


H' 


Zljv (ci ; • ■ • CiV ) B JSf ^ • • ■ 5 ) -^iV (cKni; • ■ • 


^(^)n 


where we have defined 


Biv(a„,,...,a„„) = —-- / d^a; det [Ai(a;j - . 

TT.-i Ai'(-a„,) Jwm l<^.J<N 


n^i Ai'(-a„J Jwk 
Similarly, the propagator in the denominator goes like 

f d^xG^°\x,T\e,0)c^T-^yMFN 

JW^f 


n< 


Ajv(£i, ■ • ■ Eiv) 


where Fn is a constant with no simple expression. 

For reflecting boundary conditions we have instead the following expression for the numerator 
f d^xG^j^\x,T\e,0) ='^ f d^xF^\x)Fn\e)e~^^^^'^ 

Jwn „ 

= ’ - - - ^ {2xr/^ 

-^v . ^ Pm ■ • • Pnjv 

L 


d^x- 


^det^__ [Ai((2A)^/^a;j - /3„J] 

7 ■ • ■ 

(dn-i ■ ■ ■ 


/V - ^ Pm • • • Pn^ 

-- -- ^ - det [AUzj—pm)] 

= ^-^A^ixl,..., x%){2X)^^ Y: 

- Pni • • • Pnjv 


L 


1 


ni=i Ai(-/3ni) Jw, 


[ , det [Ai( 2 :j - )] 

Jw„ l<t,3<N 


(34) 


(35) 


(36) 


(37) 


D On the normalisation constants 

Here we derive the expressions for the two normalisation constants for the case of reunions for abs orbi ng and 
reflecting boundary conditions. For both cases we recall the well-known result of Selberg’s integral |10| : 


1 

TV! 


roo ^oo ^ N 1 

- •••/ \A{xi,...,XNp’'Y{xYe-~'^"dxi= Y{ 

■ 7o do 


F[a + jc)F{{j + l)c) 


r{c) 


For the constant Bn we do the change of variables gf/2 = yi so that qidqi = dyt or dqi 


_ d-Vi 


(38) 

so we can write 


Bn = 


POO PO 

Jo Jo 


dqi ■ ■ ■ dqNC 




Aiv(^l5 • ■ • 5 to) 


/ py \ dy POO POO 

= (§) ...^ dyr...dy.ne--y-/^ 


Y[y^ 


4iir(l/i,--- ,1/iv) (39) 


9iAr(2Ar+l) JV-1 , , 9iiV(2iV+l) N-1 , 

n c (I + d CO +1) - n C(2 + ;)C (I + 


2=0 


2=0 
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Similarly for the case of reflecting boundary conditions we have 


Fiv = ( - 
, TT 


= A^!- 


N 


1 2jv(iv-i) 


2N/2 
N{2N-1) jV-1 


roo PC 

Jo Jo 


T-/2^-yi 


tN 


j=0 


dyi---dyNA%{yi,--- ,yN)Yly, 

N(2N-1) jV-1 


+ j]r{j + i) = 


n ^(9 +j)^o'+2) 

i=o 


(40) 


E Derivatives of F{x,j) 

Starting from the following expression for F{x,'y)'. 

V3 


F{x,^) = — ‘^u^^{x,^)e Q, !^,ri(a;, 7 )V u(a;, 7 ) 

X^-R VD J / 


27^ 

27*2 


and using the following property 

dU{a,b, x) 


dx 


= {a — b + x)U{a,b,x) — U{a — l,b,x) . 


(41) 


(42) 


we notice that we can write the expression 


for some set of coefficients still to be determined. The expression (431 is certainly correct for n = 1 and 

n = 2. Let us them assume is holds for any n and performe one more derivative with respect to n: 


9*"+i *"+2y3)( 




(n) 


(-7-3n + 6 n(*, 7));7 - -l,-,u{x,'y) -6u{x,^)U - -£,-,«(*, 7 ) 


/ n 


6 


(44) 


But using the property (421, we can write 




1 


6 


u{x,^)U V - -,m(*,7) ] = [r-^-o+'^]U[--£, -,u(a;, 7 ) ] - U [ - - £- 1, -,«(*, 7 ) . (45) 


6 


6 


Gathering results we find 


d F{x,'y) _ \/Z xj^^'^ix 


^=0 


(n) 


(-n + 2^);7 - -l,-,u{x,i) ] +2U !,-,«(*, 7 ) 


1 


6 


(46) 


On the other hand, we want to write this result as (431 for n —>■ n + 1. This implies to rewrite the sum as: 

n r / -1 A \ /l 


e^o 


(—n + 2£)U [ R — £, -,«(*, 7 ) ) +2(7 - — £ — I, -,«(*, 7 ) 


6 


= pC^"H-n + 2£)U (^^-£,^,u{x,j)^ +2pcl^^U (^^-£-l,^,u{x,^)^ 
= E (-n + 27)[/ Q - 7, ^, u(*, 7 )) + 2 ct\u (^ ^, u(*, 7 )) 

= Xic:'"+'>[/Q-£,^,u(*, 7 )) . 


( 47 ) 
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This results in the following set of recurrence relations for the set of coefhcents 

=C'f^(-n + 2£) + 2C'i:>i, (48) 

^(n+1) _ ry^in) 

^n + 1 ~ 5 


with the initial condition = 1. By checking explicitly the value of some of these coeffcients, and with the 
help of the Sloane databas^ we arrive at the solution: 


d"’ = 


2 "-2'^(n-£)!(2£-n)! ’ 


' = 0 ,..., n . 


(49) 


A Similar analysis can be perform ny doing derivatives with respect to 7 . Starting from (431 one notices that 

a'=+’*F(a;, 7 ) 


d^kdx” 




Ti K / ^ . 

i=o s=o 


(50) 


for some set of coefficients (f^)}- Performing one more derivative with respect to 7 allows us to arrive at 
the following set of recurrence relations for those, viz. 

= + + D^o\e), 

Di’‘+^\l) = -(^^ + k + 3£ + 3s^Di'^\£)-3Di%{£), s = l,...,fc, (^1) 

with the initial condition Hq°^(£) = 1. As we are only interested in the orders k — 0 and fe = 2, we do not need 
a general solution. This yields 

Di°'>{£) = l, 

Di^'>{£) = -^il + 2£) , d[^'>{£) = -3 (52) 

= ^(1 + 21?)(5 + 6 ^), L»d(f() = 3(7 + 6 ^), Di^'>{£)=9. 


F Moments 

We first recall the following identity 


From here we write 


dM+ 1 ) 


POO 

I 1 u —xs 

/ dss^e = 

Jo 




(l + Ai) 


POO 

/ F{x) 

Jo 

-1 POO POO 

— -— / dxF{x) / dss^e 

r{p + 1) Jo Jo 

-1 POO POO 

ni^)Jo "^'Jo 

F 7 T 1 )/ 

1 r 

r{fj, + 1 ) Jq 




2/3 


(53) 


(54) 


2 


At https://oeis.org 
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But 






Jo 

3 r 

2 Jo 

3 -3, 

2 ' 

3 _3 
^ 2 “ 


dyyi('^+®)+l/2g-«!/ 


= £„-f(M+B+l) 


[ di/y5('^+'®)+i/2g-a 

Jo 


(55) 


^'‘+-®+^'r0(At + £; + i)) 


Thus denoting = 1 + /i we finally have 

M_, 


= 3 -|(^+B) -^ (|(^ + -^)) 

2 r(:.) 


(56) 
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